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Abstract 

Star formation rates (SFRs) obtained via extinction corrected Ha are compared to dense gas as traced 
by 12 CO (J = 3-2) emission at the centers of nearby galaxies, observed with the ASTE telescope. It 
is found that, although many of the observed positions are dusty and therefore heavily absorbed at Ha, 
the SFR shows a striking correlation with dense gas in the form of the Schmidt law with an index 1.0. 
The correlation is also compared between gas traced by 12 CO (J = 1-0) and application of Ha extinction 
correction. We find that dense gas produces a far better correlation with SFR in view of surface density 
values. 
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1. Introduction 

The present knowledge of star formation on galactic 
scales in relation with its precursor gas is generally ex- 
pressed by the Schmidt law (Schmidt 1959), 

SFR(xp N (1) 

where SFR is the star formation rate, p the gas den- 
sity, and N the Schmidt law index, expressing the effi- 
ciency of star formation from gas. Often written also in 
terms of surface averaged quantities, equation (1) relates 
two physical values SFR and p which are generally spa- 
tially decoupled when observed locally, and connected in 
a spatially averaged sense. The connection between the 
two values also has a time averaged nature, namely the 
formation timescale of massive stars. Therefore, in order 
to obtain valid physical suggestions from the Schmidt law, 
we must derive these two values based on measurements 
which express conditions that are spatially and temporally 
connected as much as observations allow. 

Previous observations of molecular gas have been con- 
ducted mainly in the 12 CO (J = 1-0) line, tracing 
cold gas which constitutes the bulk of galactic molecu- 
lar clouds. However, Kohno et al. (1999) have shown that 
a tracer of denser gas, HCN, shows a better spatial cor- 
relation with star forming regions. Observations of dense 
gas tracers such as HCN (Gao and Solomon 2004) and 
higher transition CO lines like 12 CO (J = 2-1) (Braine 
et al. 1993; Boker et al. 2003) and 12 CO (J = 3-2) 
(Mauersberger et al. 1999; Yao et al. 2003; Narayanan et 
al. 2005) are becoming accessible, and it is of interest how 
these tracers of dense show up in terms of the Schmidt 



law. 

The calibration of SFR becomes an issue under these 
circumstances, where virtually all observational inquiry 
of the Schmidt law using dense gas has resorted to FIR as 
the SFR calibrator. This in part owes to the fact that FIR 
data from the IRAS satellite is abundant, and that many 
of the sample galaxies observed in dense gas tracers (HCN 
and higher CO transitions) were selected according to FIR 
luminosity. However, we must bear in mind that the SFR 
must be calibrated using massive stars. The SFR derived 
from FIR luminosity can only trace star formation over 
~ 10 8 years, because of its contamination from extended 
dust heated by the interstellar radiation field (Kennicutt 
1998), which can amount to a significant fraction of the 
FIR luminosity (e.g., Hirashita et al. 2003). 

This results in overestimation of the SFR, and also be- 
comes a bottleneck in improving the temporal connection 
of gas and star formation. Another practical shorthand 
is that the angular resolution of IRAS data is commonly 
several arcminutes; too large to infer physics from the cor- 
relation between dense gas tracer data, which are typically 
tens of arcseconds in angular resolution. The most reliable 
massive SFR tracer to date is the Ha luminosity tracing 
star formation over several 10 6 years, whose correlation 
with dense gas, surprisingly has not been checked. The 
main reason is that Ha is weak in dusty galaxies such as 
those observed in dense gas tracers, due to the extinction 
within the galaxies. It is important therefore not only to 
acquire the Ha luminosity of these galaxies, but to accu- 
rately estimate the amount of extinction. 

The main objective and result of this Letter is to exam- 
ine the correlation between 12 CO (J = 3-2) tracing warm 
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dense gas (typically ~ 30K), and extinction corrected Ha 
luminosity tracing accurately the SFR, both in surface 
averaged densities. 

2. Observation 

Observation of the 12 CO (J = 3-2) at 345GHz was 
conducted using the Atacama Submillimetcr Telescope 
Experiment (ASTE) (Ezawa et al. 2004; Kohno 2005), 
a 10m single dish located in the Atacama desert of al- 
titude 4800m in Pampa La Bola, Chile. Observations 
were remotely made from an ASTE operation room of 
the National Astronomical Observatory of Japan (NAOJ) 
at Mitaka, Japan, using a network observation system N- 
COSMOS3 developed by NAOJ (Kamazaki et al. 2005). 

The sample was selected to be able to compare the dense 
gas quantity and its relation with star formation, in both 
normal and starbursting galaxies. The galaxies were se- 
lected so that most them have readily accessible and ex- 
tinction correctable Ha data. Another limitation for the 
samples was that their velocity width had to be under 
350 km s _1 , chosen so that the emission fits within the 
bandwidth 445 km s _1 of the backend, allowing for base- 
line subtraction. Table 1 lists the observed samples. The 
galaxies were observed only at their central position, and 
the resolution FWHM of ASTE (22 ") corresponds to a 
linear size of ~3 kpc for a typical sample distance of 30 
Mpc. 

The observations were conducted in August 2005, in fair 
weather conditions, using a double sideband cooled SIS 
mixer. Calibration was done using the standard chop- 
per wheel method. Backend was a 1024 channel digi- 
tal spectrometer with 512 MHz bandwidth and frequency 
resolution of 0.5 MHz, corresponding to a velocity reso- 
lution of 0.43km s _1 . The resulting velocity bandwidth 
is 445km s _1 . Typical system temperatures at 345 GHz 
ranged from 180K to 300K. Pointing was checked every 
several hours using Uranus or Mars, and was found to be 
accurate to ~ 1". 

Obtained data were reduced with NEWSTAR, an AIPS 
based software, used commonly at Nobeyama Radio 
Observatory (NRO). After flagging bad spectra, first to 
second order baselines were subtracted, then smoothed to 
a velocity resolution of typically 15 to 20 km s _1 . T^ was 
then converted to T m b via T m b = T^/0.6, were 0.6 is the 
main beam efficiency of the ASTE telescope. The resul- 
tant integrated intensity of the galaxies are listed in table 
1. 

3. Results 

The integrated intensity Ico(^ = 3-2)= /T m bdw can 
be converted to surface gas density using the conversion 
factor Xqo such that S oc Acolco, assuming that 12 CO 
( J = 3 2) uniquely traces dense gas. However, gas tem- 
perature can contribute to 12 CO (J = 3-2) emission and 
will introduce complications. In order to circumvent this 
difficulty, we use Ico (J = 3-2) hereafter, regardless of 
what it implies physically. Errors for the values listed in 
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table 1 were calculated using 

SIco = o-^AVcoSV [Kkms- 1 ] (2) 

where a is the r.m.s. noise in T m b, AVco the full line 
width, and SV the velocity resolution (15 or 20 km s _1 ). 
The obtained spectra are presented in figure 1. 

All galaxies were detected except for NGC 520, but 
we attribute this to NGC 520's line width of 500km s" 1 
Solomon et al. (1992). 

3.1. SFR 

The SFR was derived for the sample galaxies where pos- 
sible, using narrow band Ha imaging data from Young et 
al. (1996). Flux within the ASTE beam was calculated 
using task "apphot" on the FITS images using IRAF soft- 
ware. The SFR surface density Esfr was then calculated 
using the formulation by Kennicutt (1998), and corrected 
for inclination by a factor cosi. 

Ssfr (M q pc~ 2 yr" 1 ) 

= 7.9xHr 42 ^f^cosz (ergs s-^c" 2 ) (3) 
o 

where S is the projected area of the observing beam. The 
SFR derived in this way was then corrected for internal 
extinction into Xgp^ when possible, by dust using E(B — 
V) magnitudes derived from Ha/H(3 ratio given in Ho 
et al. (1997). Combined with equation 3, we obtain the 
extinction corrected SFR SgpR. 

In order to compensate for the lack of data, we com- 
piled 12 CO (J = 3-2) from other sources. Table 2 lists 
these sources, along with the SFR derived as above where 
possible. 

3.2. Schmidt Law 

Figure 2 shows the obtained Schmidt law between ex- 
tinction corrected Ha and 12 CO (J = 3-2). For com- 
parison, we also show the relation between 12 CO (J = 
1-0) and Ha. 12 CO (J = 1-0) data are taken from sim- 
ilar resolution (16" or 22") surveys by Komugi et al. (in 
preparation), Braine et al. (1993), and Nishiyama and 
Nakai (2001). Apparently, a combination of extinction 
corrected SFR and dense gas gives a better correlation. 
Table 3 gives the Schmidt law index N and correlation 
coefficient. In all cases, the index N of the Schmidt law is 
found to be effectively 1, consistent with previous studies 
which use total luminosity in comparing the two values 
representing gas and SFR (Gao and Solomon 2004; Yao 
et al. 2003; Boker et al. 2003), even though we use sur- 
face densities which are more indicative of the intrinsic 
properties of the samples. 

3.3. Inclination Correction 

An effect we must consider when using surface aver- 
aged values is the inclination of the target galaxy, where 
the surface density is derived by multiplying the luminos- 
ity by cosi. However, inclination is generally defined from 
its global morphology, and hence is not always indicative 
of the disk in its central kpc, as in this case. The cen- 
tral region may be thick or clumpy, so that corrections for 
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Col.(l):Galaxy name. Col.(2)(3):Coordinates from NED. Many were observed also at optically defined coordinates 
from Dressel & Condon (1976), used in the 12 CO (J = 1-0) survey by Komugi et al. (in preparation) , marked * . For 
NGC 157 and NGC 925, coordinates from both NED and Dressel & Condon (1976) match. Col.(4)(5)(6):Morphology, 
distance, and inclination in degrees, as listed in Young et al. (1995). Col. (7) Observed integrated intensity of 12 CO (J 
= 3-2), converted to main beam temperature units. Col. (4) (5): SFR and extinction corrected SFR, respectively, as 
explained in text. 



Table 2. Sample Galaxies. 



Table 3. Least Squares fit 
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Col.(l): Galaxy name. Col. (2): References. 1 refers 
to Mauersberger et al. (1999), 2 refers to Narayanan et 
al. (2005). Both were observed at the HHT, with an an- 
gular resolution of 22", same as ASTE. Col. (3): Integrated 
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Col.(l): Tracer used for dense gas. Col. (2): "No" for no 
Ha extinction correction, "Yes" for correction applied as 
explained in text. Col. (3) (4): Schmidt law index TV and 
correlation coefficient, from a least squares fitting. For all 
gas tracers, the lower row is for no inclination correction. 



inclination are not necessary in the first place. To see if 
a factor of cosi is justified, we plot the inclination cor- 
rected gas densities and SFR vs. cosi in the left hand 
side of figure 3. Apparently, the values decrease with in- 
creasing inclination (i.e., as they become edge-on), a trend 
which should not be seen if the correction is correct, be- 
cause inclination of galaxies is a meaningless parameter 
regarding its intrinsic properties, and the average face-on 
surface density should be independent of inclination. The 
decreasing trend is roughly the same for 12 CO (J = 1-0), 
12 CO (J = 3-2) and SFR, which suggests that this trend 
is not due to opacity of the emission lines, but more of a 
geometric quality, i.e., that the structure inside the 22" 
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Fig. 1. 12 CO (J = 3—2) spectra of the galaxy cen- 
ters observed at ASTE. Spectra with an asterisk arc 
those observed at positions as explained in table 1. 



Fig. 2. The obtained Schmidt law for the sample galaxies, 
using combinations of gas ( 12 CO (J = 1-0) and 12 CO (J 
= 3—2)), and SFR with and without internal extinction cor- 
rection. The lower right hand panel is best correlated, with 
N = 1.0 shown as the best fit line. 

beam is more clumpy than disky. The resulting Schmidt 
law will be biased towards N = 1, because the gas and SFR 
densities decrease equally along a slope of 1. The right 
hand side of figure 3, is the same for values uncorrected 
for inclination. We do not see any trend with inclination, 
which is physically more plausible. 

From these views, we conclude that a correction for in- 
clination of cosi is not completely justified; figure 4 show 
the Schmidt law for all values uncorrected for inclination. 
The variance of the residuals for a best fit is given in table 

3, and in this case also we see that the use of 12 CO (J = 
3-2) and extinction corrected Ha improves the fit greatly. 

4. Discussion 

We have shown for the first time that 12 CO (J = 3-2) 
has a striking correlation with Ha derived SFR. Assuming 
that 12 CO (J = 3-2) indeed traces dense gas, this implies 
that Ha may be a valid SFR tracer even in dense and 
dusty regions. Dense gas 12 CO (J = 3-2) also correlates 
better with recent star formation compared to 12 CO (J = 
1-0). 

Interpreting this result in a qualitative manner is easy. 
Assuming that star formation occurs where gas density 
exceeds a certain value, we can expect that 12 CO (J = 
3-2) is more spatially and temporally connected to star 
formation compared to 12 CO (J = 1-0). By using Ha as 
a SFR tracer, the spatial connection (resolution) and tem- 
poral connection (traces SF over 10 6 years) are even more 
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Fig. 4. Same as figure 2, but uncorrected for inclination. 
1 CO (J = 3—2) with corrected SFR (correlation coefficient 
0.78, N=) still shows a better correlation compared to 12 CO 
(J = 1—0) (correlation coefficient 0.41). 
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Fig. 3. Inclination versus surface gas density and SFR den- 
sity, corrected for (left) and uncorrected for inclination 
(right). Notice that figures on the left show decreasing density 
with inclination, whereas figures on the right show no trend. 

improved. This improvement should show up in terms of 
the Schmidt law. 

Using Ha. we have circumvented the need to use global 
gas and SFR in expressing the Schmidt law, which is 
known to introduce a size effect. Larger galaxies tend to 
be more luminous in any wavelength: therefore produc- 
ing a correlation between gas and SFR without any other 
physical reasons (see Stark et al. 1986). The resolution 
of Ha imaging allowed us to meaningfully derive surface 
density values, whereas using the IRAS data may force us 
to compare spatially decoupled gas and SFR within the 
large beam. 

The Schmidt law index N, however, should be treated 
with care. The correction for extinction can be applied 
only to detectable emission, and not the intrinsic total 
Ha emission. This is not unique to recombination lines, 
but to all star formation tracers where extinction plays 
a role. In that light, The observed Schmidt law index 
is always an underestimate, where in our case the effect 
may be strong because our study focuses on the central 
regions of dusty galaxies. This is in concord with the 
literature, where the widely accepted value of 1.4 (e.g., 
Kcnnicutt 1998) is considerably higher than our derived 
value of ~ 1.0 
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